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Abstract 

We discuss the collective excitations in a spatially inhomogeneous (cylindrically symmetric) Bose- 
Einstein condensation (BEC) at low temperature (T <C T\). The main result is the dispersion 
relation for a (first) sound wave that is obtained by describing the perturbation as a Lagrangian 
coordinate. The dispersion curve is in good agreement with the Bogoliubov phonon spectrum 
to = ck, where k = kg = m/r, the wave number and c = [AirapTi 2 ] 1 ! 2 /M , the speed of first 
sound. Based on Bohm's quantum theory, a spontaneously broken gauge symmetry in a quantum 
fluid is discussed in terms of the quantum fluctuation-dissipation, from which it is shown that the 
symmetry breaking takes place at the free surface of BEC in an external field. 
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I. INTRODUCTION 



The Bose-Einstein condensation (BEC) is a remarkable quantum phenomena occurring 
in a macroscopic Bose system l|. London has proposed that the A transition between He 
I and He II is a result of the same process which causes the condensation of an ideal Bose 

Q 

gas at low temperatures |2J. In 1941, Landau studied the collective excitation spectrum in 

He II based on his phenomenological theory of two-fluid model and has shown that the 

excitation spectrum u(k) of a wave number k should rise linearly with slope c [i.e., u = ck 

with the speed of (first) sound c = 238m / sec] for small k, pass through a maximum, drop to 

a local minimum at some value k , and then rise again for k > ko, which has been confirmed 

fl 

in great detail by Henshaw and Woods |4|. For small k, the excitations are called phonons 
(quantized sound wave). 

In the absence of a satisfactory microscopic theory, there has been considerable devel- 
opment of phenomenological theories following Landau's two fluid model. Of the various 
attempts to account for the A transition, the best known is that of Feynman j^j. Based on 
the exact quantum partition function as an integral over particle trajectories by his path- 
integral approach to quantum mechanics, Feynman studied the two-fluid model of He II in a 
strongly interacting 4 He gas and pointed out that London's view on BEC is essentially cor- 
rect. The point in his argument is that, in a liquid like quantum system, the symmetry of a 
Bose system plays far more important a role than that of strong inter-atomic forces which do 
not prevent these particles move freely in the system; yet the pair-interaction which brings 
about A-transition also ensures that He II possess the collective excitations (phonons and 
rotons) |5[. It is also the essential point of Bogoliubov's theory of superfluidity which has 
laid the basis for much of our theoretical understanding about superfluidity of He II 00. 

Penrose and Onsager further extended the concept of Bose condensation to a strongly 
interacting Bose gas of He II based on the properties of a ground state wave function derived 
from the assumption that there is no long-range order interaction. They proposed that one 
could essentially define the Bose condensed state of a system such as He II as a state in 
which the reduced density matrix of the system can be factorized in a certain special way, 
now known as the off-diagonal long, range order (ODLRO) after Yang jjj. 

Although the previous theories 0, 0, [j| have been very successful in explaining the collec- 
tive phenomena in He II, the dispersion relation for a longitudinal (first) sound wave u = ck 



2 



is given only in the case of a spatially homogeneous system for mathematical convenience [t| • 
Bogoliubov j(| has used an ingenious mathematical method based directly on quantum field 
theory for his calculation of the excitation spectrum in a weakly interacting Bose system. 

The formal extension of the Bogoliubov method to a finite inhomogeneous system is 
difficult, since one can quantize a scalar field of quasi-particles only in a Hilbert space 
Furthermore, it is impossible to make use of the elegant mathematical technique developed 
in his quantum field theory approach, because the theory of collective excitations in a finite 
space problem reduces to an initial-boundary value problem for which one must find a correct 
boundary condition in terms of a dynamical variable. It is, therefore, important to use an 
appropriate mathematical technique to develop a theory of collective excitations in BEC in 
a trap. 

Recently, however, the study of a new BEC in a trap has helped to revive interest in 
the study of collective excitations, but it has faced an almost insurmountable challenge to 
derive a correct dispersion relation for collective excitations in BEC in a trap [ljj. The 
problem with the boundary conditions appearing in the dynamical equations represents the 
major stumbling block in obtaining a proper dispersion relation for phonon spectrum at low 
temperature (T <C T A ) 

With proper boundary conditions a theory must be capable of computing the excitation 
spectrum for phonons J(| , and yet explains in a natural way why a symmetry braking takes 
place and how the sound wave dissipates by the fluctuation-dissipation process at the free 
surface. These are uniquely new phenomena in a Bose system that is confined by an external 
potential. Such a theory is the subject of this paper. 

In this paper we consider an imperfect Bose gas of N identical particles in a trap. Between 
each pair of any particles in the system, there is a hard sphere repulsion of range a that 
balances the external potential which determines the size of the condensation in the trap. 
This pair-interaction is the basic, and almost the only, assumption to describe the weakly 
interacting Bose particles in a trap. Moreover, this assumption is reasonable to give results 
which are satisfactory in a dilute Bose gas in trap ^j] ; the results are also consistent with the 
theory of superfluidity [6] . The main contributions of the present paper are a semi-classical 
calculation of the excitation spectrum for phonons and a study of the symmetry breaking 
along with the fluctuation-dissipation process of sound waves to conserve the energy in an 
isolated system. The study of collective excitations is essential for a complete understanding 
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of certain properties of a Bose condensed system, in particular, the longitudinal collective 
excitations (first sound) at temperature near absolute zero. This provides the course of 
experiments that confirm BEC in a trap Q]. 

Because of a finite, inhomogeneous density of BEC, the collective excitation spectrum 
depends on the spatial density profile and the boundary conditions. In this paper we shall 
study the existence and excitation energy spectrum of a cylindrical (first) sound wave in a 
cylindrically symmetric inhomogeneous Bose condensate (CSIBC) which is an extension of 
the previous study [18]. The present analysis, in principle, can be extended to an ellipsoidal 
Bose condensate in a trap. In practice, however, the mathematical analysis of the collective 
excitations in an ellipsoidal Bose-Einstein condensation in prolate spheroidal coordinates is 
mathematically more challenging due to the complications of the ellipsoidal wave propaga- 
tions, and the well-known special cases may also play an important role to check the general 



case 
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2l|. 



To present this work as a self-contained paper, we first review the concept of ODLRO in 
which Penrose and Onsager suggested that the reduced density matrix of a Bose condensed 
state can be factorized, 

p(r, r') — ip'(r)tp(r) + j(\r — r'\. (1) 

where 7 — > as \r — r'\ — > 00. The single particle wave function if)(r) represents the 
condensed state in ODLRO and is viewed as a function of macroscopic dynamical variables. 
The single particle state function ip(r) is also interpreted as the mean value of a quantum 
particle field [lol ]. 

With the hard sphere approximation, one can show that the mean field satisfies the 
nonlinear Schrodinger equation (Gross-Pitaevskii), 

r)ih ft 2 

iH lK = ~2M V ^ + [V{XU + 9ll ^> (2) 
where V se if = gi\^p\ 2 = Airh 2 a/M\ijj\ 2 and a is the s-wave scattering length. It is a self- 
consistent Hartree equation for the Bose condensed wave function ll| . In our model of 
an imperfect Bose gas, the wave function describes a near perfect gas modified as little as 
possible by the presence of the pair-interactions. Basic to the theory of superfluidity is the 
idea of Bogoliubov [f| that this pair-interaction which brings about the Bose condensed state 
also ensures that the system possess the longitudinal collective excitations (first sound) in 
BEC. We must therefore use the effective potential that includes the self-interaction term in 



Schrodinger equation, which might be thought of as a generalization of Bohm's theory ^| 
to the theory of superfluidity p, since the nonlinear term ^i|V'| 2 is invariant under a U(l) 



group transformation 
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II. BASIC EQUATIONS 

In order to derive the density profile for CSIBC, we write the equations for the ensemble 
average energy in the usual quantum theory Q|: 

U = J V> f ("|^V 2 + V(x) ext + |m 2 ) v<fa, (3) 

£ave = S (^ iv ^ 2 + v ^ x)ext ^ 2 + f 1 ^' 4 ) dx > S ^ dx = N ' (4) 

For a BEC in a trap, the ground state density can be obtained in terms of an external 
potential and the chemical potential by minimizing the energy functional Eq. (jlj) for a fixed 
number of particles in the system with the condition p = (Penrose-Onsager criterion for 
BEC): 

p(x) = mx)\ 2 = -IL-fr - V(x) ext ), (5) 
Ann a 

where ft is a Lagrangian multiplier and is the chemical potential. It should be noted that 
the ground-state wave function Eq. (J3j) has a nodal surface at which the density becomes 
zero and is the boundary of CSIBC. 

As we have emphasized above, the most serious question in regard to a possible application 
of Bogoliubbov's quantum field theory technique (or an extended version, the Bogoliubov- 
de Gennes equation) to a finite Bose system is that the theory remains valid only in an 
infinite uniform system or a system in a box with periodic boundary conditions but not a 
finite inhomogeneous system of BEC in a trap [t| ■ 

Both for this reason, and because it is essential to show that a symmetry breaking takes 
place at the nodal surface of CSIBC, we employ Bohm's quantum theory with emphasis 
on the ensemble of particle trajectories Q, ^j]. Briefly, Bohm writes ip(r) in the form 
: f(x, t)exp[jrS(x, t)] = p l / 2 {x, t)exp[jrS(x, t)], where S(x, t) is a phase (or an action) 
We then rewrite Eq. (j2J) to obtain, 

dp „ , VS, , , 



Mx, t) -- 
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dS (VS) 2 ,„ . h 2 ,V'p 1 (Vp) 2 . 

where V(x) = V(x) ext + gp. Here the potential V(x) includes V se if = gp- The last term 
Ueqmp of Eq. (j6b|) is the effective quantum-mechanical potential (EQMP) which plays a 
crucial role in the discussion of a symmetry breaking in a finite spatially inhomogeneous 
Bose system. 

In Eq. (|B|L Bohm suggested that one may regard Eq. (jHajl as the conservation of current 
if one identifies p as the probability density and v = V S(x) /M . S(x,t) the solution of 
Eq. (j6b|) is the action. However, in the limit h — > 0, S(x,t) is a solution of the Hamilton- 
Jacobi equation and becomes a phase of the wave function ip(x,t) 14,^^. In general, the 
solution of the quantum Hamilton- Jacobi equation Eq. (|6bJ) defines an ensemble of possible 
particle trajectories, which can be obtained in principle from the Hamilton- Jacobi function, 
S(x), by integrating the velocity, v(x) = V S(x) / 'M . The equation for S implies, however, 
that the particle is acted on, not only by a potential V(x) but also by the effective quantum 
mechanical potential U eqmp . Moreover Eq. (}6b"j) suggests us how we choose stable particle 
trajectories about which we may linearize the equations of motion. For if, as was done in 
the previous study of collective excitations in He II j^O], a system is described as an infinite 
uniform medium, U eqmp becomes zero. On the other hand, in a finite system the effective 
quantum mechanical potential diverges at the nodal surface, 

n ft 2 r vy i{v P )\ _^vv 

eqmp 4M [ p 2 p 2 M / ' 1 } 



where p(r) is zero at the nodal surface. 
As emphasized by Bohm 



l^ . a particle experiences the force from U eqmp and fluctuates 
with its momentum p = VS* and energy near the surface with the degree of divergence 
Ueqmp = Mu^po/ D 2 as D = [p — (l/2)Mu 2 r 2 ] — > near the surface. Hence U eqmp breaks 
up the phase coherence of a sound wave in the surface layer, the thickness of which is in the 
order of mean free path. 

III. FIRST-ORDER EQUATIONS 

Now we come to the question of why we must find an alternate approach such as Bohm's 
quantum theory to the quantum field theory method by which Bogoliubov has already ob- 
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tained a correct energy spectrum of collective excitations u = ck with c = [An aph 2 ] 1 / 2 / M 
in the long wavelength limit. The reason is obvious: one can quantize the scalar field of 
quasi-particles only in a Hilbert space jflQ], but not in a finite inhomogeneous system such 
as CSIBC. By the same reason, an extended version of the Bogoliubov theory with the 
mean field, the Bogoliubov-de Gennes equation, is not also applicable to a finite inhomo- 
geneous system and it would lead to erroneous results. This observation provides a deeper 
appreciation of Feynman's path-integral approach in his atomic theory of the two-fluid model 
although his analysis was confined to a uniform system |5| . The basic approach in both Feyn- 
man's two-fluid model and our present approach is that a Bose system with self-consistent 
interactions can be studied by the ensemble of particle trajectories to which Bohm's quantum 
theory is more useful for describing a system with a well-defined boundary. More impor- 
tantly, by Bohm's interpretation of quantum mechanics, the solution of the Hamilton- Jacobi 
equation can be separated from the solution of the quantum Hamilton- Jacobi equation in the 
classical limit (i.e., h — > 0), and the action of the quantum Hamilton- Jacobi equation is used 
only for the measurement process to reconcile with the uncertainty principle by statistical 
ensemble of particle trajectories. This unique feature allows one to apply perturbations to 
the solution of the Hamilton- Jacobi equation, and at the same time the effective quantum 
mechanical potential is useful to describing the fluctuation-dissipation process to explain the 
spontaneously broken local gauge symmetry in the system. 

As in the previous work 18|, the perturbation to a particle trajectory is treated as a 
Lagrangian coordinate to the semi-classical solution. This leads a set of the linearized 
dynamical equations from which we derive a second-order partial differential equation in 
terms of a displacement vector. We solve this differential equation for the solution of a 
problem in CSIBC with correct initial-boundary conditions. 

Since the present calculation of collective excitation spectrum in CSIBC is new, and is 
of some interest in itself, we present it here in full detail. We now introduce a symmetry 
breaking perturbation to the particle trajectories in CSIBC in a Lagrangian coordinate 
23l 124 125| ; it is defined as x = x + £(x ,t), where £(x ,t) is a function of the unperturbed 
position of a particle in the condensate, and remains attached to the particle as it moves. 

In addition to the particle displacements, we introduce phase coherence that relates the 
Bose condensed wave function (mean field) in ODLRO to a many-body ground-state wave 
function, ^■'VS(x ,t) = X]« & ' ^^cm^cm; where S 0ti (x 0;i ,t) is the phase of a single Bose 
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particle in the system. This is a necessary and sufficient condition for the phase coherence. 
It also gives a simple interpretation of the mean field and quasi-particles (coherent excited 
states) in BEC. It is therefore obvious that a single particle excitation does not occur in our 
study as in Feynman's work [5|. Moreover, the collective excitations we study in this paper 
are phase coherent sound waves. 

The first step in developing this alternate approach is to associate with a Bose particle 
having precise stable particle trajectories which are a function of position and momentum. 
In this connection, it is worthwhile to note that the use of the Hamilton- Jacobi equation in 
solving for the motion of a particle is only matter of convenience and that, in principle, we 
can always solve it directly by the conventional dynamical equations with the correct initial- 
boundary conditions for a finite space problem. In order to study the modified excitations 
in a trap, we find it convenient to write down the basic dynamical equations. 

This is most simply done by writing v(x , t) = VS(x , t)/M, then Eqs. © and © yield 
the following dynamical equations: 

the equation of motion for a single particle 

d 

M(—v + v-Vv) = -Vp, (8) 

the equation of continuity, 

l p + V -(pv) = 0, (9) 

where p is henceforth interpreted as the number density, 
the equation of state, 

p(x, t) = H\ oc \p{x, t)\ + V exk , (10) 

n 

where pi oc [p(x, t)} is the chemical potential in the local density approximation [26]. These 
dynamical equations allow us to show that the symmetry breaking takes place at the nodal 
surface and yield a correct modified phonon spectrum in CSIBC. 

A little algebra with Eqs. (jH)), © an d (fTUj) gives the following linearized equations: 



S(x,t) =S(x ,t)+$- V S(x ,t) (11a) 

v(x,t) = ^t (lib) 

p(x, t) = p(x ) - V • \p(x Q )$] = p(x ) + 5p (11c) 

d 

fl{x,t) = Ploc(x ) + V cxt - 7^(/iloc)[V • (M)]> ( lld ) 
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where p{xq) = in BEC and Vo denotes the partial derivative with respect to Xq with 
V -> V - V £ ■ V . Eq. (fTTcl) was derived by substituting Eq. (Illbjl to Eq. ©, then 
integrating over time. 

These full equations of motion describe the ensemble of stable particle trajectories from 
which we may the study of the perturbed particle trajectories as we have done previously 
for a spherically symmetric condensate Q|. The perturbed particle orbits are obtained by 
the same method as Weinberg's analysis of electron orbits 2q . 

It may also be worth of noticing, though less obviously, that Feynman's atomic displace- 
ments with the back-flow in his intuitive description of phonons in 4 He p are mathematically 
equivalent to our approach with the Lagrangian displacements and the phase coherence in 
ODLRO; we arrive at the almost same dispersion relation for a sound (longitudinal) wave. 
Hence our analysis of collective excitations in BEC is essentially equivalent to that of atomic 
theory of the two-fluid model of Feynman since both approaches study the ensemble of par- 
ticle trajectories. 

In order to make a contact with a possible future experiment, we make the following 
assumptions: the condensate is a long and thin cylinder, i.e., L> b, where L is the length of 
the condensate, b the radius. We therefore assume kb <C 1, where k is the wave number along 
the z-axis. We may relax the last assumption, but with considerable amount of additional 
algebra. We further assume, for simplicity, a cylindrically symmetric condensate and A = 
oj° z /(ju\ <C 1 which is consistent with L>6, and also fc6 < 1. These assumptions reduce the 
problem effectively to a two-dimensional problem in the cylindrical polar coordinates. 

If we take the linear momentum po = in BEC, the equation of motion gives 

/iioco(PoOo), t) + Kxt = constant = /xi oc o[po(0)] = Po, (12) 

which shows that the peak density at the center of a trap plays a pivotal role - an essential 
point, since the speed of (first) sound c = [Att ap(O)^ 2 ] 1 / 2 / 'M must be finite at the center and 
dissipates to zero giving rise to a surface energy at the nodal surface as we shall see. 

Now we are ready to derive the first-order equation of motion in terms of £, which de- 
scribes the collective excitations with initial-boundary conditions. To show this we linearize 
Eq. © using Eq. (JTHJ), Eqs. jnj) and Eq. and the result is 

|U = JjVa - u>* U ■ V)r + (r ■ V£)] - ^ 2 [or + \r 2 Va]. (13) 
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Here we have taken V ext (r) = McUq r 2 /2, dpi oc /dp = Aitfi 2 a/M |l2|, a = V ■ £, and co is 
the radial trap frequency, and have also dropped the subscript in r . It is also understood 
henceforth that r = r±, V = Vj_, and | = (£ n t;g)e ime . 

The entire discussion of collective excitations is based on the first-order equation Eq. (TEH) 
which must be solved by initial-boundary conditions [18]. We proceed further in two stages. 
First, we discuss the surface waves on the free surface of CSIBC in a trap and show that 
the free surface of the BEC in equilibrium behaves like a classical fluid. Secondly, we solve 
the same equation, but remove the condition of incompressibility to derive the dispersion 
relation for a compressional wave, an ordinary sound wave (phonons) which shows that 
the presence of collective excitations with excitation energies that are linear u = ck with 
respect to both the speed of first sound c = [An aph 2 ] 1 ^ 2 / M , and the phonon wave number 
k = kg = m/r, consistent with the Bogoliubov theory p. 

IV. DISPERSION RELATION FOR SURFACE WAVES IN CSIBC 

In this section we carry out the first part of the analysis for the collective excitations in 
CSIBC. The excitation of surface waves on a free surface can be initiated by perturbations 
of the external trapping potential just like gravity waves on the surface of a fluid in a 
gravitational field. 

For the surface waves, we may impose the following two boundary conditions on £ at the 
free surface of the condensate: one of which is incompressibility of a fluid at the free surface 
and the other of which is the condition of irrotational motion at the free surface 0, 2(3 ] . 
Since we intend to relax the first condition later in the discussion of a compressional wave 
which is our main objective in this paper, no extensive discussion is necessary. In any case, 
they are mathematically o = V • £ = and uj = V x £ = 0. These two conditions allow 
one to rewrite the first-order equation, Eq. (JT3Jl as 

|^ = -wg*-w8(r-V)£ (14) 

Also, the above two conditions imply that there exists a scalar function x that satisfies 
the Laplace equation, V 2 x = and £ s = — V%. Here the subscript s stands for the surface 
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waves. The general solution for \ is given in the cylindrical coordinates by 



x (r, t) = [Q™(t)r m + Q™(t)r- {m ~>} e imd (15) 

where we may set Q™{t) = for a solid cylindrical condensate. 

It is a simple algebra, taking the gradient on Eq. (JT5Jl and substituting it into Eq. (JHJ), 
to obtain the time-dependent equation for Q™, 

^Q™ + mu 2 Q™ = 0, (16) 
from which one can write down the dispersion relation at once 

U 2 S urf = mUJ l ( 17 ) 

The dispersion relation for the surface waves is a function of the radial trapping frequency 
(external force) ujq and the mode number m. There are two aspects to the dispersion 
relation. First, the dispersion relation is independent of the the pair-interaction potential 
(i.e., internal dynamics) of a trapped Bose gas and of the radius of the condensate. A priori, 
the condition of irrotational flow is assumed to derive the dispersion relation. Since this 



condition is valid for a surface wave both in a classical fluid and in a superfluid j20L |22|, 
the dispersion relation Eq. (J 17)) is a manifestation of broken symmetry in a Bose system. 
Moreover, it is independent of % in spite of the quantum ground state density given by Eq. (J3J) 
which is a function of h. Second, the surface waves are driven solely by the perturbations of 
an external trapping force on the surface of BEC just like a gravity wave under the action 
of the force of gravity in a gravitational field j^]. 

Unlike the spherical wave , the dispersion relation Eq. ()17|) has an interesting geomet- 
rical interpretation: for m — 0, it corresponds the sausage mode for which the fluid column 
makes a radial oscillation, but this mode is absent in CSIBC on account of to = 0; for m — 1, 
the hose mode, the column makes a garden hose like motion without changing the shape 
of the cross-section. For higher mode m ^> 1, kg — m/b corresponds a short waves along 
the circumference of the column with the radius of b. Moreover, it is fairly simple a task 
to observe the surface waves by applying the perturbation along the waist line in a long 
ellipsoidal condensate. The experimental data [27[ agree well with Eq. (fT7|) . 

A few final remarks will close this section. First, it is important to note that the dispersion 
relation for the surface waves in BEC [Eq. (fTTj) ] also implies that Osborne's experimental 
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observation on the contour of a rotating He II 



2 a ] remains indeed correct, a point which 



appears to be quite contrary to what one expects from Landau's two-fluid model j[j and 
which has been extensively studied j^J [ji]] since Osborne's first observation on the surface 
contour of rotating 4 He |28| in 1950. In short, the free surface of a superfluid behaves 
like a classical fluid under the external force; this is a peculiarly universal behavior of a 
superfluid. It is also evident that Feynman's vortex line model j^] cannot explain the 
Osborne's observation as emphasized by Meservey [3^. Second, Rayfield and Reif 2 
observed that, apart from the quantization of circulation, a vortex moves like a classical 
fluid with an ion probe at core, which implies a break-down of superfluidity at a vortex 
core This peculiar universal classical fluid-like behavior of the free surface of CSIBC is 
essentially in agreement with the breakdown of superfluidity at a nodal surface of a vortex 
core. This observation also resolves the recent controversy over the Magnus force (3^. The 
classical Magnus force employed by Vinen JsfJ in his analysis of vortex quantization is indeed 
the correct one which gave the value k — $ v s ■ dl, where k = h/M = 0.997 x 10 _3 cm 2 /sec 
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V. DISPERSION RELATION FOR (FIRST) SOUND WAVES IN CSIBC 

One problem that requires the introduction of the collective excitations for its solution 
is that of the superfluidity of CSIBC. Our aim here is to understand the (first) sound 
wave propagation and its modified excitation spectrum in CSIBC. This calculation gives the 
excitation spectrum almost identical to that of Bogoliubov's phonon spectrum, u = ck with 
c = [47r aph 2 ] 1/2 /M and k = m/r. 

It may be useful to recall Feynman's theory of phonons at this point. In a series of papers 
Q], Feynman has laid out an elegantly simple theory of collective excitations (phonons and 
rotons) in 4 He based on his path-integral approach to quantum mechanics based on the exact 
quantum partition function. Feynman has shown that, at temperature sufficiently below 
the A point, the ground state wave function has a positive amplitude for any configuration, 
since the ground state wave function is symmetric and thus has no nodal surface. A density 
fluctuation involving a large number of particles (or long wavelengths) creates a back-flow 
by the conservation of particle density; the compressed density in one part of the system is 
then left with the rarefaction adjacently. Feynman goes on to argue that, with the back-flow 
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and the conservation of the momentum, the atomic displacements do not lead to a single 
particle excitation, but to extreme low-energy excitations of compressional waves with no 
nodes (phonons), since a nodal surface breaks the symmetry property of a Bose system as we 
have shown in Sec II of this paper. In spite of a simple picture of sound waves in Feynman's 
analysis just like Landau's two-fluid model jjjj], the speed of (first) sound is not given as a 
function of the pair- interaction potential which he presented at the outset, but is shown to 
be identical with the velocity of sound defined by the usual macroscopic considerations [^. 

In contrast, the above description of our method implies that, without any reference to 
a classical fluid, one should be able to explicitly calculate the modified excitation spectrum 
as a function of the speed of sound and the wave number with the condition of irrotational 
motion of a superfluid. We now define the compressibility of the fluid as a = V • £ in 
Eq. (jllcj) . It is a mathematical description of small-amplitude density waves (phonons) 
through the equation of continuity and the phase-coherence, and is essentially equivalent 
to Feynman's picture of phonons in superfluid He II. The energy spectrum of collective 
excitations in the phonon regime can be obtained from Eq. ()13j) by taking the divergence on 
both sides. With the vector identity V • [(r • V)£] = r • Vcr + a, where V = Vj_ and r = rj_, 
and the condition of irrotational motion V x £ = [3j, it is a straightforward algebra to 
obtain the first-order equation, 

J^<r(r, t) = \4 {a 2 - r 2 )VV - 3 uj 2 r |-a - 4 uj 2 a, (18) 

where a 2 = 8tc an (0)h 2 / (M lu ) 2 and V 2 = V^. 

Writing a(r, t) = S(t)W(r)e , we obtain the variable separated equations: 

^S(t) + X n S(t) = 0. (19) 

(« 2 - r 2 ) [l±( r ±W n (r)) - ^ W n (r)} 
r ar ar r z 

-6r—W n (r) 

+2(X n /uj 2 -A)W n (r)=0, (20) 

where A„ is a constant of separation. 

The dispersion relation is determined by the eigenvalues of the radial equation, Eq. (j2Uj). 
In the following we show the eigenvalues are a function of the speed of (first) sound and the 
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radial trap frequency by transforming Eq. (|20|) to the Sturm-Liouville problem and thereby 
obtaining the eigenvalues in terms of the complete set of ortho-normal functions, 

An = [ b r(a 2 - r 2 f [(^-W n ) 2 + % W 2 (r)}dr (21) 
Jo dr r 

and with 

f r{a 2 -r 2 ) 2 W m {r)W n {r)dr = 5 m , n . (22) 
Jo 

Since X n = 2(\ n /uo 2 - 4) and b = [8vr ap (0)h 2 ] 1/2 / (u M), Eq. d2U> shows that the eigen- 
values are indeed a function of the speed of (first) sound at the center of the condensate 
(i.e., b = [v2 / 'uo\c ctr and c ctr = [Anap^h 2 ] 1 / 2 / M) for all values of the mode number m. To 
derive Eq. (|2ip. we have tacitly assumed a point source at the center of CSIBC and thus 
cr(r) = at both r = and the free surface r = b. Eqs. ({21} and (|2*2*|) show that a sound wave 
travels from the center, where it has the peak velocity, to the free surface since the speed of 
the (first) sound wave varies as the density of CSIBC is not uniform. The eigenvalues are 
also a function of the mode number. The reason for this is of course that the mode number 
m is a constant of the motion. 

In principle the eigenvalues can be obtained from Eq. 1)2 ip. but the integral representa- 
tion is hopelessly complicated to evaluate the eigenvalues. Nevertheless, it shows that the 
eigenvalues are positive definite and they are indeed a function of the speed of (first) sound 
c ctr for all values of the mode number. 

It is helpful to look at Eq. (j2"U|) in complex plane and to consider Eq. (j2"Uj) as a linear 
eigenvalue equation. This can be done easily by mathematical transformation; it also serves 
a very useful and important purpose to see how the sound wave travels in CSIBC. If we define 
the complex variable W n (r) = r ±{ ^Z n (r) and with x = r 2 /a 2 , then Eq. (|27]jl transforms to 
a well-known differential equation, 

d 2 d 
x(l-x) ^Z n + [c-(a + b + l)x] —Z n -abZ n = 0. (23) 
dx l dx 

This is the Gauss's hypergeometric equation 37J with c = ±m + 1, a + b = ±m + 3, and 
ab = (1/4) [A n =F 6m]. 

The advantage of this transformation is to calculate the eigenvalues by means of a simple 
numerical method. Moreover, there are none of the end-point singularity problems associated 
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with W m (r) that one encounters in the integral representation Eq. (|21|). And the various 
numerical methods [38] are available to compute the eigenvalues. 

One further point must be made with regard to Eq. (|23|) . The hypergeometric equation 
Eq. (J23j) has regular singular points at x — 0, x — 1 and x = oo. Its solution is the 
hypergeometric function which is analytic in the complex plane with a cut from 1 to oo 
along the real axis. The branch point x = 1 corresponds to the location of the free surface of 
the fluid column of CSIBC. Parenthetically, we also note that the analyticity of the complex 
variable reflects the fact that a sound wave can travel only inside of the free surface as it 
must. 

Since the present work critically depends on the numerical accuracy, we may pause to 
examine to what extent one can take the numerical solutions to the eigenvalue problem. 
One question that arises is whether or not it is possible to see the physical significance of 
the cut in the complex plane. To better understand the analyticity, it should be emphasized 
that cr = V • £ has been interpreted as an amplitude of a compressional wave (i.e., collective 
mode) in CSIBC in which a cylindrically symmetric sound wave travels from the the center 
toward the free surface. The presence of the cut implies that there cannot be a sound wave 
beyond the branch point. 

The collective solution exists in a domain in which o ^ 0. Such a solution must be 
analytic in the complex plane. To check this we must be sure a solution to Eq. (J2*3|) is 
compatible with Eq. (}2T| . It is physically obvious that a sound wave can travel in the 
domain [0, 1] in which the solution to Eq. ()23j) is analytic. The branch point reflects the 
nodal surface at which the sound wave must dissipate to conserve the total energy of an 
isolated system giving rise a surface energy. Thus our numerical solutions in the domain 
[0, 1] will always satisfy this physical condition; it is a unique feature of finite systems. 

Here we employ a simple numerical method j^] to evaluate the smallest eigenvalues cor- 
responding to the low-energy phonon excitations in the domain [0, 1]. This simple numerical 
method should, however, be quite satisfactory provided the solutions are numerically stable 
near the branch point 38]. 

Returning to Eq. (j!9|) and taking S(t) = e tuJt , we obtain the dispersion relation as 

u ph = ±(X n ) 1/2 = ±[2X s + 4} 1 / 2 uj . (24) 
where X s is the smallest eigenvalues from Eq. (j2HJ). The ratio uj/ujq with respect to the mode 
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number m is plotted in Fig.l. 

The dispersion relation Eq. (j23| (Fig. 1.) shows that for large m, uj/ujq is linear with 
respect to the mode number m (or wave number kg = m/r). It is almost exactly the same 

nn 

form of the usual dispersion relation u = ck of a (first) sound wave [3j, |6J with kg = m/r. 
This has the following simple interpretation: for large values of the mode number m, a wave, 
whose wavelength is shorter than the size of CSIBC but longer than atomic dimensions, sees 
the medium as if it is an infinite medium. Thus we recover in an essential manner the 
Bogoliubov dispersion relation u = ck (f| at long wavelengths. This assures us that in 
fact the present analysis describes the collective excitations correctly at temperature near 
absolutely zero in CSIBC. For small values of the mode number, the dispersion curve is 
slightly parabolic with respect to the mode number due to geometrical effects-finite size 
effects. For m = mode, lo/ujq — 2.8153 is a unique value in a finite space problem. It 
corresponds to a uniform perturbation in polar angle in CSIBC. 

As in the case of a spherically symmetric nonuniform condensate, the whole problem of 
collective excitations in CSIBC is reduced to the question of how one measures the speed of 
first sound c c t r = [47ran(0)fr 2 ] 1//2 /M (or the peak density) for an outgoing cylindrical sound 
wave at various radial points in CSIBC. This requires a precise measurement of the peak 
density, a prerequisite for an experimental confirmation of BEC in a trap. 

There is one more point to be made about observing the phonon spectrum in a spatially 
inhomogeneous condensate. To understand the difficulty, let us take the simplest case m = 
mode. Unlike the surface waves, it is important to perturb the medium at the center to 
initiate a (first) sound wave uniformly in polar angle. How, then, if we were to confirm 
BEC in CSIBC by studying the radial propagation of the sound wave, might we try to to 
initiate short waves at the center and to measure the speed of sound wave at r < b. As far 
as we can see, this is a difficult task; it is, nonetheless, essential to see if we can initiate 
the sound waves at the center and place a probe at r < b, inside the nodal surface, since 
the the sound waves dissipate near the nodal surface as discussed above. The question 
of whether a systematic perturbation consistent with the present theory is experimentally 
feasible remains to be seen. Nevertheless, it would be necessary to have a precise dispersion 
curve on the collective excitations in CSIBC similar to that of Henshaw and Woods in He II 
And it is still the only way to quantitatively establish the realization of BEC in a trap. 
There is one last point to be made about the break-down of superfluidity in BEC at 
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the free surface. To understand this, let us consider what has happened to the sound 
wave propagating toward the free surface. As it approaches the free surface, the particle 
trajectories rapidly fluctuate due the effective quantum mechanical potential U P „ m r, which 
diverges as Mto 2 fj, /D 2 as D = [fj, — (l/2)Mulr 2 \ — > near the free surface |l4|. This 
fluctuation breaks up the phase coherence of the sound wave which dissipates at the free 
surface giving rise a surface energy. 

The important mechanism for the symmetry breaking is that the sound wave dissipates 
by the interaction with the particles in the surface layer which is ensured by the fact that 
the initial state and the final state of CSIBC always have the same energy. Therefore, the 
p_on of the fluctuation of particles due to V m and the dissipation of sound waves 
at the surface layer can be understood in terms of Kubo's fluctuation-dissipation theory |40J. 
One of our main objectives of this paper is to understand how the symmetry of the ground 
state wave function breaks down at the free surface, i.e., a nodal surface. It is natural 
to identify this underlying basic mechanism for the symmetry breaking as a spontaneously 
broken symmetry at the free surface which accompanies phonons as Nambu-Goldstone bosons 
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VI. CONCLUSION AND DISCUSSION 



I should like to close with one final remark on the symmetry breaking in an isolated BEC 
in trap. By describing the perturbation as a Lagrangian coordinate, our analysis on the 
dynamical equations is a accurate quantum mechanical treatment for a many-body problem 
just as that of the quantum field theory method by Bogoliubov for a uniform system 
provided a correct quantum ground state function Eq. (jSJ) is given. Not only does the 
analysis correctly predicts the low-energy excitation spectrum for phonons, consistent with 
the Bogoliubov spectrum but also shows the correct low-energy spectrum for the surface 
waves in a finite inhomogeneous system confined by the external potential. 

The basis of our reasoning for the symmetry breaking at the nodal surface follows from 
the observation that, in spite of the quantum ground state given by Eq. the dispersion 
relation for the surface wave Eq. (|T7jl is independent of h but not for the phonon spectrum 
of a superfluid Eq. (J24|) . This is a definite proof that the gauge symmetry of a Bose system is 
broken at the free surface. In the phonon regime the Bogoliubov dispersion relation uj = ck 



17 



holds in CSIBC with the small geometrical corrections. One also notices from Fig.l that 
there is a striking similarity between the dispersion curves of the spherical and the cylindrical 
condensates, although the angular momentum i and the mode number m have an entirely 



different meaning 25|. The salient point in our results is that the modified energy spectrum 
for the phonons is a function of the speed of first sound, c ctr = [Anan{Q)h 2 ] l l 2 / M (or the 
s-wave scattering length). As Bogoliubov emphasized in his work on a uniform Bose gas at 
zero temperature (|, the weakly repulsive pair- interaction potential indeed plays a crucial 
role for the structure of the ground state and the superfluidity 
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FIG. 1: The ratio uj/lo = A 1 / 2 = [2A + 4] 1 / 2 is plotted clS el function of the mode number m. It 
shows how the energy spectrum of phonons varies with the mode number. 
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